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On  the  dissipation  associated  with  equilibrium 
shocks  in  finite  elasticity^ 


by 


James  K.  Knowles 

Division  of  Engineering  and  Applied  Science 
California  Institute  of  Technology 
Pasadena,  California 


Abstract 


'Equilibrium  fields  with  discontinuous  displacement  gradients  can 
occur  in  finite  elasticity  for  certain  materials.  The  presence  of  such 
^equilibrium  shocks^  affects  the  energy  balance  in  the  elastostatic  field, 
and  the  present  paper  is  concerned  with  a notion  of  dissipation  associated 
with  this  energy  balance.  A dissipation  inequality  is  proposed  for  three- 
dimensional  equilibrium  shocks  for  both  compressible  and  incompressible 
materials.  The  consequences  of  this  inequality  are  studied  for  weak  shocks 
in  plane  strain  for  compressible  materials  and  for  shocks  of  arbitrary 
strength  in  anti-plane  strain  for  a class  of  incompressible  materials.  A 
thermodynamic  argument  for  the  dissipation  inequality  is  also  given. 

1 . Introduction 

An  example  discussed  in  [1]  shows  that  it  is  possible  for  the  diffe- 
rential equations  governing  finite  equilibrium  deformations  of  an  elastic 
solid  — for  an  otherwise  "reasonable"  material  — to  lose  their  ellipticity 


The  results  communicated  in  this  paper  were  obtained  in  the  course  of  an 
investigation  supported  by  Contract  N00014-75-C-01  96  between  the 
California  Institute  of  Technology  and  the  Office  of  Naval  Research. 
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in  the  presence  of  sufficiently  severe  strains.  For  a homogeneous,  isotropic, 
compressible  elastic  solid,  explicit  conditions  on  the  principal  stretches 
and  the  strain  energy  which  are  necessary  and  sufficient  for  the  ellipticity 
of  the  field  equations  of  plane  elastostatics  were  given  in  [2]. 

One  of  the  consequences  of  such  a loss  of  ellipticity  is  the  possible 
occurrence  of  elastostatic  fields  in  which  the  displacements  may  fail  — in 
certain  portions  of  the  body,  at  least  — to  have  continuous  derivatives  of 
even  the  first  order.  A detailed  investigation  of  the  local  structure  of  such 
fields  in  plane  finite  elastostatics  is  described  in  [3],  where  it  is  pointed 
out  that  they  may  be  relevant  to  the  study  of  the  "Liiders  bands”  commonly 
observed  in  ductile  materials*.  Indeed,  deformations  with  abruptly  changing 
gradients  arise  in  a wide  variety  of  solids  in  connection  with  what  has  been 
termed  "localized  shear"^.  Boundary  value  problems  in  the  equilibrium 
theory  of  finite  elasticity  in  which  locally  severe  deformations  — near  a 
cavity,  for  example  — give  rise  to  a loss  of  ellipticity  and  an  associated 
reduction  in  the  smoothness  of  the  stress  and  displacement  fields  would 
seem  to  be  of  special  interest  in  the  study  of  fracture. 

The  mathematical  description  of  phenomena  of  the  kind  described 
above  has  certain  features  in  common  with  the  theory  of  stationary  tran- 
sonic  flows  in  inviscid  gas  dynamics^.  One  of  the  conspicuous  aspects  of 
such  flows  is  the  occurrence  of  shocks  — surfaces  across  which  the  fluid 

^See  Chapter  18  of  [4]. 

^For  an  extensive  discussion  of  localized  shear,  see  the  article  by  Rice  [5], 

^In  fact,  a precise  analogy  exists  between  the  boundary  value  problems 
for  gas  flows  past  an  obstacle  on  the  one  hand,  and  finite  anti-plane  strain 
(for  a certain  class  of  incompressible  elastic  materials)  of  a medium  con- 
taining a crack  or  cavity  on  the  other  hand.  See  § 7 of  [6], 
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velocity,  pressure,  density  and  entropy  suffer  jump  discontinuities.  In  the 

setting  of  finite  elastostatics , the  analogous  surfaces  are  those  which  carry 

jump  discontinuities  tn  the  first  derivatives  of  the  displacement  vector; 

they  will  be  called  equilibrium  shocks  in  the  present  paper. 

In  gas  dynamics,  one  of  the  essential  restrictions  on  the  flow  across 

a shock  — whether  moving  or  stationary  — is  that  imposed  by  the  requirement 

that  the  entropy  of  a particle  shall  increase  as  the  particle  crosses  a shock^. 

2 

This  condition  leads  to  major  qualitative  conclusions  concerning  the  flow  , 
and  it  also  expresses  — in  an  idealized  way,  to  be  sure  — the  dissipative 
character  of  the  process  of  shock  formation. 

In  this  extensive  study  of  weak  solutions  of  hyperbolic  systems  of 

3 

conservation  laws.  Lax  has  emphasized  the  mathematical  role  of  "entropy 
conditions"  of  various  kinds  in  connection  with  the  issue  of  singling  out  the 
physically  meaningful  weak  solution  to  the  initial  value  problem  from  among 
the  many  such  solutions  nominally  admitted  by  the  differential  equations 
and  initial  conditions  themselves.  The  many  beautiful  results  pertaining 
to  this  question  which  are  discribed  in  [8]  - [11]  appear  to  be  most  appro- 
priate for  dynamical  problems  — the  propagation  of  shock  waves,  for  example. 
They  do  not  seem  readily  interpretable  in  the  context  of  finite  elastic  equili- 
brium, in  which  time  is  not  one  of  the  independent  variables  and  the  initial 
value  problem  is  of  no  apparent  significance.  The  objective  of  the  present 
paper  is  the  derivation  and  discussion  of  a dissipation  condition  which  is 

^For  an  extensive  discussion  of  shocks  in  compressible  flow,  see  [7]. 

2 

Gas  dynamical  shocks  are  always  compressive,  for  example. 

3 

See  [8]  - [11],  where  references  to  related  work  may  be  found. 
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analogous  to  the  "increasing  entropy"  requirement  lor  steady  gas  flows  and 
which  is  thought  to  be  appropriate  for  equilibrium  shocks  in  finite  elasticity. 
A detailed  application  of  the  result  derived  here  may  be  found  in  [3j. 

The  field  equations  of  finite  elastostatic s for  both  compressible  and 
incompressible  materials  are  given  in  the  following  section.  The  class  of 
weak  solutions  to  be  considered  is  descriced  in  Section  3,  which  also  in- 
cludes the  jumpconditions  to  be  satisfied  across  an  equilibrium  shock.  The 
dissipation  condition  for  a quasi- static  time-dependent  family  of  equilibrium 
shocks  is  derived  and  discussed  in  Section  4.  The  implications  of  the 
dissipation  condition  for  weak  shocks  in  plane  deformations  of  compressible 
materials  are  examined  in  Section  5.  Section  6 is  devoted  to  equilibrium 
shocks  and  the  dissipation  inequality  in  '■.nti-plane  strain  for  a special  class 
of  incompressible  materials.  A thermodynamic  argument  in  support  of  the 
dissipation  requirement  is  given  in  the  final  section. 


2.  Equations  of  elastostatics 

Let  R be  the  interior  of  a three-dimensional  region  occupied  by  an 
elastic  body  in  its  undeformed  state.  In  a deformation  of  the  body,  R is 
mapped  invertibly  onto  a domain  R , so  that  a particle  with  position  vector 
X in  R is  carried  to  the  point  in  R whose  position  vector  is 


y = y(x)  = X + u(x)  , xeR  . (2.  1) 

For  the  moment,  the  displacement  vector  ju  is  assumed  to  be  twice  con- 
tinuously differentiable  on  R. 

Let  x^,  y^  and  u^  be  the  components  in  a fixed  rectangular  cartesian 
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coordinate  system  of  x,  y and  ja,  respectively.  The  tensors^  F,  G and  C 
defined  by 


F=7^=J^  + 7u  , G=FF^  , C-F^F 


(<i.  2) 


where  7 is  the  gradient  operator  and  1 the  identity  tensor,  have  components 

.3 


F. . = y.  . = 6 . + u.  . , G..  = F , F.,  , C..  = F,  .F  . ; 

ij  ^i,j  ij  i,j  ij  ik  jk  ij  ki  kj 


(2.  3)- 


F is  the  deformation  gradient  tensor,  G and  C the  left-  and  right-  defor- 
mation tensors,  respectively.  The  Jacobian  of  the  mapping  is  assumed 
positive; 

J = det  F >0  on  rt  . (2.  4) 


If  p and  Py  are  the  mass  densities^  in  R and  rt  balance  of  mass 


requires 


pJ  = p. 


on  R 


(2.  5) 


If  Ojj  are  the  components  of  the  nominal  (or  Piola)  stress  tensor 
the  local  condition  of  force  equilibrium  in  the  absence  of  body  forces 


^ The  same  symbol,  e.  g.  F,  will  be  used  for  a tensor  and  its  matrix  of 
components  in  the  given  coordinate  system. 

^The  superscript  T stands  for  transposition. 

Latin  subscripts  have  the  range  1,2,3  and  repeated  subscripts  are  summed 
over  this  range.  A subscript  preceded  by  a comma  indicates  differentiation 
with  respect  to  the  corresponding  x-  coordinate. 

4 

Pq  is  taken  to  be  constant. 

The  o.j  represent  forces  per  unit  undeformed  area. 


6 


may  be  expressed  as 


a.  . = 0 

ij.j 


on  rt  , 


(^.  6) 


provided  that  O is  continuously  differentiable  on  . The  components  T 
~ ij 

of  the  true  (or  Cauchy)  stress  tensor^  T are  related  to  C . via  F as 

~ ij 

follows: 


T.  . — O.,  F., 

ij  J ik  jk 


If  is  continuously  differentiable  on  R , then  (2.6)  is  equivalent  to 


(^.  7) 


3T. 


-LL  - 


avj 


(2.  8) 


Moment  equilibrium  leads  to  = ^ji'  ^ general  not  symmetric. 

For  a compressible  homogeneous  elastic  material  with  a strain  energy 
density  W = W'(F)  per  unit  undeformed  volume,  the  relation  between  stress 
and  deformation  may  be  taken  in  the  form 


o =iwtF) 

ij  ar.. 

■'  ij 


(2.  9) 


or,  equivalently  because  of  (2.7), 


T - i 


The  represent  forces  per  unit  deformed  area. 


(2.  10) 
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The  displacement  equations  of  equilibrium  for  a compressible  material 
then  follow  from  (2.  9),  (2.  6)  and  the  first  of  (2.  3)  as 


9 awcF) 

9x.  ar. . 

J iJ 


ijkf 


on  , 


(2.  11) 


in  which 


a^W(F) 

aF..aF, . 
IJ  kf 


(2.  12)^ 


For  an  isotropic  compressible  elastic  material,  W depends  on  F 
only  through  the  invariants  and  J : W = W(Ij,l2,J),  where 

Ij=TrG  . - 1/2  [(TrG)^  - Tr  (G^)]  , (2.13) 

and  Tr  denotes  the  trace.  In  this  case  (2.  9)  specializes  to 


a. . 

U 


ik^kj 


+ 


Mf- 

aj  ^ji  ' 


(2.  14) 


where  F.  is  the  i.j  - element  of  the  inverse  F ^ of  F, 

IJ  ~ ~ 

given  by  the  last  of  (2.  3). 

For  an  incompressible  homogeneous  elastic  material, 
deformations  which  preserve  volume  locally  are  admissible, 
(2,  10)  are  replaced  respectively  by 


and  C. . is 
ij 

only  those 
and  (2.  9), 


^W(F)  is  assumed  to  be  infinitely  differentiable  for  every  nonsingular 
tensor  F. 
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1 


» 


a. . 
ij 


9W(F) 

ar.  ~ 

ij 


(^.  15) 


T.  . 

IJ 


9W(F) 


aF. 


ik 


U.  16) 


where  is  the  Kronecker  delta  and  p is  an  arbitrary  hydrostatic  pressure, 

assumed  to  be  continuously  differentiable  on  R . The  displacement  equations 
of  equilibrium  for  an  incompressible  material  are  now 

together  with  the  constraint  expressing  incompressibility: 

J = detF  = 1 on  R . (2.  18) 


The  functions  (2.17)  are  again  given  by  (2.12);  (2.17)  follows  from 

(2.6),  (2.  15)  and  the  fact  that  F'^  . = 0 if  detF  = 1.  The  system  (2.  17), 

Ji » J ~ 

(2.  18)  comprises  four  equations  for  the  four  unknowns  u.  ,p. 

For  an  isotropic  incompressible  material,  W depends  on  F only 
through  the  invariants  Ij , of  (2.  13),  and  (2.  15)  specializes  to 


ij  Nolj  1 ij  ik  kj  ^ ji 


(2.  19) 


3,  Weak  solutions  and  equilibrium  shocks 

In  place  of  the  smoothness  assumptions  made  in  the  preceding  section, 
it  will  now  be  required  that  the  displacement  vector  u be  merely  continu- 
ous and  piecewise  continuously  differentiable  on  R,  while  the  nominal 
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stress  tensor  0 is  piecewise  continuous  on  R.  For  incompressible 
materials,  it  is  further  assumed  that  the  hydrostatic  pressure  p is  piece- 
wise  continuous  on  R Thus  F,  C_  (and  p,  for  incompressible  materials) 
may  suffer  jump  discontinuities  across  certain  surfaces^  in 

Enlarging  the  class  of  admissible  stress  and  displacement  fields  in 
this  way  makes  it  necessary  to  generalize  the  sense  in  which  certain  of  the 
field  equations  are  to  be  satisfied.  The  tensors  F,  G and  C whose  com- 
ponents are  defined  in  (2.  3)  are  now  piecewise  continuous  on  li,  as  is  the 
Jacobian  J of  (2.4).  Equation  (2.  5),  expressing  the  local  balance  of  mass, 
continues  to  have  meaning  under  the  present  circumstances  and  is  in  fact 
equivalent  to  the  integral  condition  expressing  mass  balance  in  the  large. 

The  local  equations  (2.6)  of  force  equilibrium,  however,  must  be  replaced 
by  the  corresponding  global  conditions  from  which  they  were  originally 
derived: 


for  every  ciosed  regular 
surface  S lying  in  rt  . 


(3.  1) 


Here  N is  the  unit  outward  normal  on  S.  The  integral  version  (3.  1)  of  the 
force  balance  requirement  remains  meaningful  for  nominal  stress  fields 
^ which  are  merely  piecewise  continuous  on  rt  . 


The  components  T_  of  true  stress  are  still  to  be  calculated  from 
and  F^j  according  to  (2.  7).  Regarded  as  functions  on  R , the  T_ 
are  piecewise  continuous.  The  symmetry  condition  T_  = continues  to 
be  equivalent  to  the  global  balance  of  moments  under  present  circumstances; 


this  imposes  via  (2.  7)  the  restriction 


* These  surfaces  are  assumed  to  be  regular  in  the  sense  of  Kellogg.  [12]. 
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Under  the  smoothness  assumptions  now  in  force,  ju  , F and  ^ are 
said  to  furnish  a weak  solution  of  the  field  equations  for  compressible 
materials  if  (2.  9),  the  first  of  (2.  3),  (2.  5),  (3.  1)  and  (3.  2)  hold.  The 
inequality  (2.  4)  is,  of  course,  still  to  be  imposed. 

For  incompressible  materials,  a weak  solution  is  furnished  by  _u, 

F,  and  p provided  the  above  equations  hold,  except  that  (2.  9)  is  to  be 
replaced  by  (2.  15). 

It  is  possible  to  enlarge  still  further  the  class  of  admissible  fields 

u,  F,  2 (and  p,  where  appropriate).  One  such  enlargement  would  permit 

2 

unbounded  displacements  and/or  stresses  at  isolated  points  or  along  curves 
in  ft.  A second  — and  considerably  more  significant  — weakening  of  the 
smoothness  restrictions  would  permit  discontinuities  in  the  displacement 
vector  u itself.  This  particular  relaxation  of  requirements  would  indeed 
be  essential  for  the  study  of  fracture.  Neither  of  these  generalizations  will 
be  considered  here. 

Suppose  now  that  the  deformation  gradient  tensor  F and  the  nominal 
stress  field  ^ associated  with  a weak  solution  (for  either  a compressible 
or  an  incompressible  material)  suffer  jump  discontinuities  across  a surface  S 
in  ft  but  are  continuously  differentiable  on  either  side  of  S in  a neighborhood 

^It  follows  from  (2.  14)  that  (3.  2)  is  automatically  satisfied  for  an  isotropic 
elastic  material.  In  fact  the  same  applies  to  all  objective  elastic  materials: 
for  these,  W depends  on  F only  through  C. 

^This  might  be  necessary,  for  example,  to  accommodate  singularities  at  the  | 

edge  of  a terminating  surface  which  carries  discontinuites  in  F or  d . | 

i 

I 


I 

i 


I 


of  a point  P on  S . Thi-  global  oquilibrivim  conditions  (3.  1)  then  furnish 

two  conclusions:  first,  by  the  familiar  argument,  u satisfies  the  local 

equilibriv-in^  conditions  (<i.  6)  on  either  side  of  S near  P.  Second,  by 

applying  (3.  1 ) to  a small  sphere  centered  at  P and  using  the  fact  that  (2.6) 

hold  on  either  side  of  S inside  the  sphere,  one  finds  that  the  jumps  in 

a.  must  be  such  that 
ij 


fo.  1'^  N.  = 0 at  P on  S . 

‘J  - J 


(3.  3) 


He  re  fa. . = o. . - a, . 

' ij  - ij  ij 


J.. 


and  a.,  stand  for  the  limiting  values  of  0.. 
»J  ij  “ ij 


on  S from  the  positive  and  negative  sides  of  S,  and  the  "positive"  side 
of  S is  the  side  into  which  N points.  Equation  (3.  3)  asserts  that  the 
nominal  traction  ^ with  components 


s. 

i 


o.  .N. 
U J 


(3.4) 


is  necessarily  continuous  across  S.  It  is  possible  to  show  that  (3.  3)  is 
equivalent  to  the  condition  that  the  true  traction  with  components 
tj  = r>.  be  continuous  across  the  deformation  image  S of  S:  here 

n is  the  unit  normal  on  S . 

Let  L be  any  unit  vector  which  is  tangent  to  S at  P.  Since  u 
is  continuous  at  P and  continuously  differentiable  on  either  side  of  S, 
near  enough  to  P it  follov/s  from  the  first  of  (2.  3)  that 

[F.  L.  = 0 at  P on  S (3.  5) 

I ij'-  J 

is  assumed  to  have  a continuously  varying  unit  normal  N near  P,  and 

the  derivatives  a..  are  presumed  to  have  finite  limiting  values  from  either 
•■J » K 

side  at  all  points  of  S which  are  sufficiently  close  to  P. 
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A surface  S in  rt  carrying  jump  discontinuities  in  F and  ^ 
satisfying  the  jump  conditions  (3.  3)  and  (3.  5)  is  called  an  equilibrium  shock. 

jjt  >(< 

To  distinguish  S in  rt  from  its  deformation  image  S in  rt  , it  is  con- 
venient  to  refer  to  S as  the  material  shock  surface  and  to  S as  the  spatial 

Jr 

shock  surface.  It  is  of  course  the  spatial  shock  S that  one  would  actually 

* 

"observe”;  nevertheless,  the  pre-image  S of  S figures  more  prominently 
in  the  present  analysis. 

To  investigate  many  of  the  questions  pertaining  to  the  local  existence 
and  local  properties  of  equilibrium  shocks,  it  is  sufficient  to  consider  the 
case  in  which  S is  a plane  through  P and  u is  such  that  F is  constant 
on  either  side  of  S.  Such  piecewise  homogeneous  deformations  and  the 
associated  equilibrium  shocks  have  been  studied  in  [3]  for  the  case  of  plane 


strain  in  compressible  materials. 


4.  Energy  considerations  and  dissipation. 

The  first  objective  of  the  oresent  section  is  the  derivation  of  a 
formula  which  reveals  the  effect  on  the  mechanical  energy  balance  in  an 
elastic  body  of  the  presence  of  an  equilibrium  shock.  I'his  formula  has- 
connections  with  a conservation  law  arising  in  tlie  elastostatic s of  homo- 
geneous materials  as  well  as  with  the  theory  of  delects  in  an  elastic  solid; 
it  will  lead  to  a notion  of  dissipation  for  equilibriun*  shocks  which  is  the 
principal  concern  of  the  present  paper. 

in  order  to  obtain  the  energy  formula,  it  is  necessary  to  consider 


It  is  proved  in  [3]  that  a necessary  condition  for  the  existence  of  an  equi- 
librium shock  in  plane  strain  of  a compressible  material  is  that  the  (plane 
version)  of  the  displacement  equations  of  equilibrium  (2.  11)  must  suffer  a 
• loss  of  strong  ellipticity  at  some  deformation.  A minor  modification  of  the 

argument  used  in  [3]  shows  that  this  result  remains  true  for  three-dimen- 
sional  equilibrium  shocks  in  compressible  materials. 

» 

1 


) 


i 
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a quasi-static  time -dependent  family  of  equilibrium  states.  Thus  let 

n (x,  t),  F(x,t),  ^(x,t)  furnish,  for  each  t in  a suitable  time  interval,  a 

weak  solution  of  the  field  equations  of  equilibrium  for  a compressible  material. 

Note  that  no  inertia  effects  are  considered  here:  the  time  t plays  only  the 

I 

role  of  a history  parameter.  Suppose  that  this  weak  solution  involves 

an  equilibrium  shock  across  the  surface  in  the  deformed  body  at 

time  t,  and  let  be  the  pre -image  in  the  undeformed  body  ft  of  S^. 

Thus  S and  S are  the  material  and  spatial  shock  surfaces,  respectively. 

3 

Let  V (x,t)  stand  for  the  velocity  vector  of  a point  on  the  moving 

surface  which  at  time  t is  located  at  x.  It  is  now  assumed  that, 

at  time  t , the  point  x on  is  such  that  V (x  , t ) neither  vanishes 

o ’ ~ o t ~ ~o  o 

o 

nor  is  tangent  to  . Thus  does  not  coincide  locally  with  a "contact 

° ° 4 

surface"  — a surface  across  which  there  is  no  flow  of  particles.  It  then 

follows  that  for  times  near  enough  to  t^,  and  for  points  x on  which, 

at  time  t , were  sufficiently  close  to  x on  , V (x,t)  will  continue 

o 

to  have  this  property. 

Now  suppose  that  £ is  a regular  subdomain  of  ft  which  contains 

X . Then  & will  be  intersected  by  for  all  times  in  a suitable  interval 
~o  ^ t 


1 


jj,  F and  a are  assumed  to  be  continuously  differentiable  with  respect 
to  t for  every  x in  ft. 

^If,  for  example,  a cylindrical  body  is  loaded  on  its  ends  to  a state  of 
uniaxial  tension,  t might  be  taken  to  be  proportional  to  the  given  applied 
force  (if  the  deformation  is  "load  controlled")  or  to  the  given  extension 
("grip  control"),  as  long  as  the  force  — or  the  extension  — is  a monotone 
function  of  real  time. 

(x,t)  is  assumed  to  be  continuous  in  (x,  t)  for  all  x on  and  all  times 
under  consideration. 

See  §52  of  [7].  Thus  (and  the’^efore  ) in  general  consist  of  different 
material  particles  at  different  times. 


I 

i 
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including  t^.  Let  denote  the  portion  of  which  lies  in  £)  . For  a 
sufficiently  small  domain  & and  for  times  t close  enough  to  t^  , one  may 

A 

choose  the  unit  normal  vector  N on  so  that  V'N  >0,  as  will  now  be 

A 

assumed.  The  positive  side  of  is  thus  the  side  into  which  V (and  there- 
fore N)  points.  Since,  to  an  observer  fixed  on  , particles  appear  to 

A 

cross  from  the  positive  to  the  negative  side,  it  is  natural  to  refer  to  the 

positive  side  as  "upstream",  the  negative  as  "downstream".  Let 

1 

be  the  open  domains  into  which  j6  is  divided  by  (see  Fig.  1). 

The  total  energy  stored  at  a time  t close  enough  to  t^  in  that  portion 
of  the  deformed  body  which  occupies  the  domain  — the  image  of  & under 
the  deformation  at  time  t — is  given  by 


U(t)  = Jw(F(x,t))dV 
D 


(4.  1) 


The  first  objective  is  to  obtain  a formula  for  the  derivative  U(t).  Clearly 


U(t)  = I W(F(x,t))dV  + jw(F(x,t))dV 


(4.  2) 


and  so  to  find  U(t)  it  is  necessary  to  calculate  the  derivatives  of  the  inte- 
grals in  (4.  2),  in  which  the  integrands  and  the  domains  of  integration  vary 
with  time.  Making  use  of  a standard  formula  for  this  purpose^  one  obtains 


U(t)  = J 


+ >J 


- 1 l^VNdA 

J Aa*  ^ 

A 


IJ 

*t 


..dV  + 1 WVNdA 
ij  J 


(4.  3) 


1 

Note  that,  while  and  depend  on  t,  & does  not, 

2 

See  [13]  , Chapter  3,  Section  6, 
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In  this  formula, 


F..  = F..(x.t)  = ^ F..(x,t)  , 
ij  ij'~  3t  ij  ~ 


(4.  4) 


and  W , W represent  the  limiting  values  of  W on  from  tile  upstream 
and  downstream  sides,  respectively. 


V = v(x,t)  = ^ u(x,t) 

'-W  ^ /V,/ 


(4.  5) 


stand  for  the  velocity  of  the  particle  which  at  time  t is  located  at  y(x,t). 


Then  from  (2.  3),  (4.  4)  and  (4.  5) 


F.  .(x,t)  = V.  .(x,  t)  . 

i,J~ 


(4.6) 


The  volume  integrals  in  (4.  3)  can  then  be  calculated  as  follows:  using  (2.  9), 
(4.6),  (2.6),  and  the  divergence  theorem. 


f 1^  F..dV  = f a.  .V.  .dV  = I a..v.v.dA 
J dF,.  ij  J ij  i,j  J ij  J I 


(4.  7) 


± U 


± ± 


where  ^ is  the  unit  outward  normal  on  the  boundary  94^  of  i> . Substituting 
from  (4.  7)  into  (4.  3)  and  observing  that,  on  that  portion  of  d&  which  coincides 

A 

with  S.  , the  normal  V = +N  , one  finds 

X 


U(t)  = [o..v.v.dA  - [(O..N.V.  + WN.V.)dA  + f (a.  .N.v.  + WN.V.)dA  , (4.8) 

' ' J ij  J 1 J ' IJ  J V i i'  J ' IJ  J i i i 

di>  s,  s. 


where  once  again  the  superscripts  + and 


on  o..  and  v.  refer  to  the 
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appropriate  limiting  values  of  these  quantities.  To  this  point  no  use  has 
been  made  of  the  jump  conditions  prevailing  across  an  equilibrium  shock. 
By  differentiating  with  respect  to  t the  equation  which  expresses  the 
equality  of  the  limiting  values  of  ^ on  S^,  one  finds  that  the  jump  in 
particle  velocity  v satisfies 


["I  “ ■ " 

in  terms  of  the  components  of  the  velocity  of  points  on  the  (material) 

shock  surface  and  the  jump  in  the  components  of  the  deformation 

gradient  tensor  F . When  (4.  9)  and  the  nominal  traction  continuity  require- 
ment (3.  3)  are  used  in  (4.  8),  it  reduces  to 


U(t)  = I s.v.dA  - 

V 1 I 

djo 


St 


(4.  10) 


where 


P.. 

ij 


W6. . 

ij 


F,  .o,  . 

ki  kj 


(4. 11) 


The  P..  are  the  components  of  the  energy -momentum  tensor 
g = - Fa,  whose  role  in  the  theory  of  defects  in  solids  has  been 

explored  extensively  by  Eshelby  [14,15,16].  Indeed,  (4.  10)  suggests  that 
[g]  51  may  be  viewed  as  a "force  per  unit  area"  associated  with  the  material 
shock  surface  . In  this  interpretation,  the  role  of  P is  precisely  that 
discussed  by  Eshelby  in  § 6 of  [16]  in  connection  with  the  "force  on  an  inter- 
face". 

In  general  the  "force  on  a defect"  has  been  shown  [14,  15,  16]  to 
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have  components 

J.  = (Sp..N.dA  , (4.  12) 

1 T ij  j 

r 

where  F is  a closed  surface  surrounding  the  defect.  If  there  are  no  defects 
inside  F,  so  that  u is  twice  continuously  differentiable  there. 


P..N.dA  = 0 ; 
J 


(4.  13)^ 


(4.  13)  represents  a conservation  law  for  compressible  homogeneous  elastic 
materials  which  was  independently  discovered  by  Rice  [17]  and  exploited  by 
him  — and  subsequently  others  — in  connection  with  crack  problems  in  elasti- 
city. The  energy -momentum  tensor  g plays  a prominent  part  in  the  dissi- 
pation condition  to  be  derived  below.  This  would  seem  to  be  related  to  the 
fact  that,  as  emphasized  by  L.ax  [ 11],  "entropy  conditions"  appropriate  for 
shocks  associated  with  the  initial  value  problem  for  a quasilinear  system  are 
connected  with  additional  conservation  laws  implied  by  the  system. 

An  alternative  and  more  convenient  form  of  the  energy  formula  (4.  10) 
is  obtained  by  noting  that,  if  L is  any  vector  tangent  to  , (4.  11)  gives 

P..N.L.  = - F,  .0,  .L.N.  , (4.14) 

ij  J I ki  kj  i J ' 


^Equation  (4.  13)  is  easily  verified  directly  with  the  help  of  the  divergence  theorem 
and  the  field  equations  of  Section  2 appropriate  to  compressible  materials. 
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so  that  by  (3.  5),  (3.  3), 

Thus  the  vector  [ P]^  N is  perpendicular  to  ; it  follows  that  (4,  10)  may 
be  written  in  the  form 

U(t)  = j£*vdA  + J [H]^V*NdA  , (4.  16) 

ai)  St 


where 


H = H(F.N)  = - P..N.N.  = - W + F .0  .N.N.  . (4.17) 

~ ~ ij  i j mi  mj  i j 

The  identity  (4.  16)  is  the  energy  formula  which  furnishes  the  basis 
for  the  notion  of  a dissipative  shock  to  be  introduced  below.  One  first  notes 
that  [H]^  =0  if  no  shock  is  present,  so  that  (4.  16)  then  reduces  to  the 
statement  that  the  rate  of  change  of  energy  stored  in  any  portion  of  the  elastic 
body  is  equal  to  the  rate  at  which  work  is  done  by  the  forces  external  to 
that  portion.  If  an  equilibrium  shock  is  present,  it  is  natural  to  require 
that  the  second  term  on  the  right  in  (4.  16)  be  nonpositive,  corresponding 
to  the  idea  that  the  shock  should  be  dissipative  in  the  sense  that  it  tend  to 
diminish  (or  at  least  not  increase)  the  stored  energy  in  the  body.  One  is 
accordingly  led  to  require  that,  at  each  instant  t , 

A 

St 


(4.  18) 


for  every  subdomain  il  of  K which  is  intersected  by  and  which  is 
small  enouizh  to  permit  N to  be  chosen  so  that  V*N  > 0 on  S = 

Suppose  (4.  18)  holds,  and  let  x be  a non-contact  point  on  , 
i.  e,  a point  at  which  V is  neither  zero  nor  tangent  to  . Since  £)  may 
be  chosen  as  the  interior  of  an  arbitrarily  small  sphere  centered  at  x , 
it  follows  that 

[H]^  ^ 0 at  all  non-contact  points  of  S^^  , (4.  19) 

provided  one  interprets  the  positive  side  of  S^  as  the  side  into  which  V 
points. 

According  to  (4.  17),  the  scalar  H in  the  dissipation  inequality 
(4.  19)  depends  not  only  on  the  particular  particle  concerned  (through  F)  but 
on  the  orientation  of  the  shock  (through  N)  as  well.  The  inequality  (4.  19) 
may  thus  be  interpreted  as  requiring  that  the  value  of  H associated  with 
a given  particle  and  a given  N shall  not  decrease  as  the  particle  crosses 
the  shock,  as  is  true  of  the  entropy  in  gas -dynamic  shocks. 

In  what  follows,  (4.  19)  will  be  required  to  hold  at  all  times  t. 

The  version  of  (4.  19)  appropriate  to  equilibrium  shocks  in  piecewise 
homogeneous  plane  deformations  of  a compressible  elastic  material  was 
derived  in  [3]  and  applied  there  to  certain  specific  strain  energy  densities. 

For  incompressible  materials,  the  argument  used  to  obtain  (4,  16) 
can  be  repeated  as  above,  except  that  (4.  7)  must  be  replaced  by 
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Here  (2.  15)  has  been  used,  as  well  as  the  relation 

F..^v.  . = dv./dy.  (4.21) 

ji  i,j  i 

in  which,  on  the  left,  is  regarded  as  a function  of  x and  t,  while  on 
the  right,  as  a function  of  ^ and  t.  The  incompressibility  condition  (2.  18), 
however,  implies  that  9v,/9y^  = 0,  so  that  the  last  integral  on  the  right  in 
(4.20)  vanishes,  reducing  (4.  20)  to  the  same  result  obtained  in  (4.  7).  It 
then  follows  that  (4.  16)  holds  also  for  incompressible  materials  with  H 
again  given  by  (4.  17).  Thus  the  inequality  (4.  19)  is  the  natural  dissipation 
requirement  to  impose  on  equilibrium  shocks  in  incompressible  materials 
as  well. 


5. 

If  the  region  R occupied  by  the  undeformed  body  is  cylindrical  with 


Compressible  materials:  weak  shocks  in  plane  strain. 


generators  parallel  to  the  - axis,  the  deformation  (2.  1)  corresponds  to 

plane  strain  if  u^  = 0 and  ^ = u^(x^  pl^tie  strain,  F,  G and 

C of  (2.  2)  are  functions  of  x,  ,x,  only  and  are  such  that  F , - F-  = 

^ 1 4 ' Qfi  30  3o 


■ ^03  ’ ^30  “ ^03  ~ ^33  <^33  = ^33  = 


Not  all  elastic  materials  can  sustain  general  states  of  plane  strain 

in  the  absence  of  body  forces.  A sufficient  condition  that  a material  shall 

have  this  property  is  that,  under  the  above  circumstances,  ^ as  determined 

by  (2.  9)  shall  be  such  that  a =0=0.  This  is  assumed  hereafter  in 

3a  o3 

the  present  section.  An  isotropic  material  always  fulfills  this  condition. 


For  plane  strain  the  displacement  equations  of  equilibrium  (2.  11) 
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in  the  compressible  case  are  readily  shown  to  reduce  to 


on  , 

o 


(5.  1) 


where  is  the  cross-section  of  the  cylinder  R which  lies  in  the  plane 

X-  = 0 , and 


c (JT)  - ^ 


(5.2) 


If  the  deformation  is  plane,  one  shows  easily  from  (2.  13)  and  (2.  4) 


that 


Ij  = 1 + I . Ij  = I + J , I = C„  = , 


(5.  3) 


30  that  for  isotropic  materials,  W = W(I,  I,  J)  = W(I,  J)  . It  can  further 

t ^ f 


be  shown^  that 


ttF«,oF,  ^ + 2jW,^(F:iF  + F^„f:^  )+  J^f:^F:^W^_  + 2W-6,^  6 
II  »p  y6  IJ\  pa  y5  ap  by/  pa  6y  JJ  I “Y 


pa-^  Y5  ap  by 

'-^F-‘ 

Pa  by 


'ay^Pb 

(5.  4) 


for  isotropic  materials;  here  subscripts  I and  J indicate  partial 
differentiation  of  W(I,J)  . 

In  plane  strain,  both  the  material  and  the  spatial  equilibrium 

shock  surfaces  are  cylindrical,  and  a material  shock  surface  intersects 

2 

the  cross-section  R in  a cui=ve  C.  The  jump  conditions  (3.3)  and 

o 


1 


See  [2],  Eq.  (1.  17). 

C is  assumed  to  have  a continuously  varying  normal. 
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1 

1 


J 


(3.  5)  become 


n: 


Np  .0 


and 


on  C 


on  C , 


(5.  5) 


(5.  6) 


where  N and  L are  the  unit  normal  and  unit  tangent  vectors  on  C,  res- 
pectively, and  the  positive  side  of  C is  the  side  into  which  N is  directed 

In  general  the  various  physical  and  geometrical  quantities  F 

± ^ 
will  vary  with  position  along  C.  Moreover,  if  the  shock 

is  regarded  as  a member  of  a quasi-static  time-dependent  family  of  equi- 
librium shocks  in  the  sense  of  the  preceding  section,  the  values  of  these 
quantities  will  also  vary  with  time.  In  the  present  section  attention  is  con- 
fined to  the  local  structure  of  the  various  fields  at  a fixed  point  on  C and 

at  a fixed  instant  of  time.  Suppose  at  such  a fixed  point  and  fixed  time, 

+ 

the  values  of  F^^p  (and  therefore  O^p)  are  regarded  as  given.  The  jump 
conditions  (5.  5)  and  (5.  6)  then  restrict  the  possible  local  values  of  F^p  , 
O^p  , and  associated  with  the  equilibrium  shock.  Since 


aw 


op  dF 


orp 


(5.7) 


and  in  view  of  the  fact  that  N and  L can  be  expressed  solely  in  terms  of 
the  local  angle  of  inclination  4*  of  L with  respect  to  the  Xj  - axis,  the 
jump  conditions  (5.  5)  and  (5.  6)  naay  be  regarded  as  four  equations  for  the 
five  unknowns  F^p  and  4>  . One  would  thus  expect  to  find  a one -parameter 
family  of  "solutions"  F^p  , 4>  of  (5.  5),  (6.  6)  for  given  F^p  . It  is  conven- 
ient to  choose  as  the  parameter  € indexing  this  family  the  relative  area- 
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change  across  the  shock,  so  tliat,  in  view  of  (2.  5), 


+ . + - 

J - J P - P 

1 ■ p 


(5.  8) 


One  now  regards  = F^p(e)  , •I'  = 4'(e)  , and  therefore  = L^(e)  and 


= N^(e)  , as  functions  of  e. 

If  in  fact  no  shock  is  present,  so  that  J J , (5.  8)  shows  that 

€ = 0.  A weak  equilibrium  shock  is  one  for  which  | F^p(e)  - | , and 

therefore  | e | , are  small  compared  to  unity.  It  is  now  assumed  that 
(5.  5)  and  (5.  6),  with  the  aid  of  (5.  7)  , determine,  for  a aiven  F , , a 

^ Qrp 

1 

smooth  family  Fj^p(e),  N^(e)  , ) near  e = 0 which  is  such  that 


"'op  • 


(5.9) 


The  objective  of  the  analysis  in  the  remainder  of  this  section  is  the  deter- 
mination of  the  weak-shock  approximation  to  the  local  value  of  the  jump 
[H]_  in  the  function  H of  (4.  17)  for  plane  strain  of  a compressible  material. 


This  necessitates  a study  of  the  corresponding  approximations^  to  F^p 


(e) 


and  N^(e). 


^Specifically  F^(e)  and  4>(e)  must  be  four  times  and  three  times  continu- 
ously differentiate,  respectively. 

^Weak  shocks  are  studied  in  detail  in  [ 3],  but  more  refined  calculations 
than  those  given  in  [3]  are  required  to  find  the  weak  shock  approximation  to 
[H]^  . 
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f 

To  begin,  one  observes  «hat  the  kinematical  jump  condition  (5.  6) 
immediately  leads  to 


(5,  10)‘ 


where  g(G)  is  an  as  yet  undetermined,  four  times  continuously  differentiable 
function  which, by  (5.9),  must  satisfy 


g„(0)  = 0 . (5.11) 

The  vector  g(c)  represents  the  jump  across  the  shock  of  the  derivative 
of  the  displacement  vector  normal  to  the  shock. 

In  order  to  analyze  the  traction  continuity  condition  (5.  5),  set 


A (G)  = 
O'  ' ' 


dwcE')  aw(^(e)) 


crp 


ar 


ap 


Np(G) 


(5. 12) 


and  observe  that,  for  all  sufficiently  small  |g  , 


A (e)  = 0 

a' 


But  (5.  12)  and  (5.2)  show  that 


a'  (e)  = 


aw(F)  aw(F(6)) 


N p(e) 


(5. 13) 


(5. 14) 


Equation  (5.  10),  which  is  not  limited  to  weak  shocks,  leads  to  the  geome- 
trical interpretation  of  the  local  deformation  associated  with  equilibrium 
shocks.  Although  this  interpretation  will  not  be  discussed  here,  it  is  ex- 
plained fully  in  [3]. 
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where  the  prime  indicates  differentiation  with  respect  to  e . From  (5.  10), 

(5.11), 


F;p(e)  = 4(€)Np(c)  , g^(€)Np(c)  , . g;(0)Np(0) 


(5. 15) 


Since  according  to  (5.  13),  '^^,(0)  - 0 , one  concludes  from  (5.  14),  (5.  15)  and 

(5.9)  that  gY(0)  must  satisfy 

where 

“„y<!J’0  = ',py6<EI'^A  (517) 

i 

f 

for  any  unit  vector  N and  nonsingular  tensor  F . It  can  be  shown  from 
(5. 8)  and  (5.  10)  that  satisfies 

= - e : (5.  18) 

differentiating  (5.  18)  with  respect  to  6 , setting  € = 0 and  using  (5.  11) 
shows  that 


g'(0)‘(F'^)^N(0)  = - 1 , (5.  19) 

so  that  (0)^  The  existence  of  a nontrivial  vector  g (0)  satisfying  (5.  16) 

requires  that  ^ 

detQ(N(0),  F)  = 0 . (5.  20) 

The  condition  (5.20)  is  necessary  and  sufficient  for  the  failure  of  ellipticity 


1 
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of  the  plane  displacement  equations  of  equilibrium's.  1 ) at  a deformation 
whose  local  gradient  is  F.  If , for  a given  F,  N(0)  satisfies  (5.  20) , then 
N(0)  is  normal  to  a characteristic  curve  of  the  system  (5.  1 ).  Foranysuch 
N(0),  one  can  show  that  (5.  16)  and  the  normalization  condition  (5.  19)  determine 
^ (0)  uniquely , provided  that  the  symmetric  tensor  Q is  not  the  null  tensor. 

It  is  unfortunately  necessary  to  compute  A^(0).  Differentiating 
(5.  14)  and  setting  € = 0 gives,  with  the  aid  of  (5. 9)  and  (5.  15), 

* N;(0)Np(0)] 

where 

a 9^w(F) 

'^apvSXpi^^^  ^ FFT  Spy^^S)  = ftF'  3F.  ’ 

Xm  ^ 0(3  y6  Xm 

Since  (5.  13)  holds  for  all  sufficiently  small  [cl,  it  follows  that  ^'(0)  = 0 

o' 

and  therefore  that  g^  (O)A^(O)  = 0.  This,  together  with  (5.  16)  and  the 
symmetry  apparent  from  (5.2),  lead  to 

^See  Section  1 of  [l] 

2 

The  case  Q * jO  is  excluded  from  consideration. 
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This  formula  will  be  used  later. 
To  compute  f H]  ^ , set 


h(e)  = [H)_^  = H(F.N(€))  - H(F(e),N(€)) 


(5.24) 


Invoking  the  special  features  of  F , £ and  N appropriate  to  plane  strain 
and  using  the  definition  (4.  17)  as  well  as  (5. 7),  one  obtains 


h(€)  = W(F(€)  ) - W(F) 


• |^^(F)N^(0F^^N^(£) 


(5.25) 


If  one  subsitutes  from  (5.  10)  for  where  it  appears  explicitly  in  (5.25), 

one  finds 


h(e)  = W(F(€))  - W(F)  - |^^(£(e))gy(e)Np(e) 


(5.  26) 


The  traction  continuity  condition  (5.  5),  together  with  (5.7),  shows  that  the 
quantity  in  brackets  in  (5.  26)  vanishes,  so  that 


h(€)  = W(E(€))  - W(F)  - |^(F(e))g  (e)N0(e) 

YP 


(5.27) 


Obviously  h(0)  = 0.  To  determine  the  small  -|e  1 approximation  to  h(€)  , I 

one  must  compute  some  derivatives  of  h at  e = 0,  Differentiating  (5.27),  | 


I 
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using  the  first  of  (5.  15)  and  referring  to  (5.  2),  one  gets 

In  view  of  (5.  II), 

h'  (0)  = 0 . 

Differentiate  (5.  28)  and  use  (5.  22)  and  the  first  of  (5.  15)  to  find 

From  (5.  11),  (5.  9)  and  (5.  15)  it  then  follows  that 

= CapY6<C>Sa(0)8Y<0)Np(0)Ng(0) 

and  hence  from  (5.  16)  that 

h"(0)  = 0 . 

Finally,  differentiate  (5.  30)  and  set  e = 0 , making  use  of  (5.  11),  (5.  9), 
(5.  15),  (5.  16)  and  the  symmetry 

h'"(0)  = 2d^pY5;^^(F)g;(0)gY(0)g((0)Np(0)N5(0)N^(0) 

+ 6‘=c»pY6<C^8i<°>8Y(0>Np(0)N'(0)  . 


(5.28) 


(5.29) 


(5.  30) 


(5.31) 


(5.  32) 


(5.  33) 
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Appealing  to  (5.  23),  one  arrives  at 


Since  h vanishes  to  second  order  at  e = 0 , one  has 


(5.34) 


h(e)  =^h'"(0)e^  + 0(€‘*)  , e-*o  , 


and  hence,  using  (5.  24),  the  weak-shock  result 


[H]'*"  = D€^  + 0(e'^)  as  e-»0  , 


(5.  35) 


(5.36) 


where,  by  (5.  33),  (5.22) 


^ 1 aV(F) 

= T7  g^(0)gY(0)g;,(0)Np(0)N5(0)N^(0) 


(5.  37) 


The  formulas  (5,  36),  (5.  37)  comprise  the  main  result  of  the  present  section. 
They  express  the  jump  [H]^  in  terms  of  the  third  gradients  of  the  strain 
energy  density  W , the  weak-shock  limiting  normal  N(0)  , and  the  weak- 
shock  approximation  g^(0)€  to  the  jump  g(€)  in  the  normal  derivative  of 
the  displacement  vector  across  the  shock;  g^O)  is  a nontrivial  solution  of 
(5.  16)  which  also  satisfies  (5.  19). 

The  fact  that,  for  weak  shocks,  the  jump  [H]^  is  of  the  third  order 
in  the  shock  strength  e is  analogous  to  the  result  in  gas  dynamics  that  the 
entropy  jump  across  a weak  shock  is  of  third  order  in  the  corresponding 


shock  strength  parameter. 


1 


1 


See  p.  143  of  [7]. 
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Consider  a quasi-static  time -dependent  family  of  plane  strain  equi- 
librium shocks  in  a given  material,  and  let  x be  the  position  vector  of  a 

point  on  the  material  shock  surface  at  a fixed  instant  t , Let  V(x,  t) 

+ + 

be  the  velocity  of  the  point  x at  the  instant  t , 51  ” 510S*^)  limiting 

value  of  51  from  the  upstream  side.  Let  N(0)  satisfy  (5.20)  and  be 

such  that  V{x»t)’N(0)  is  positive.^  Let  g^(0)  be  the  vector  determined 
uniquely  by  (5.  16),  (5.  19).  Reference  to  (5.  37)  shows  uiat  the  value  of  D 
at  the  point  x and  the  instant  t is  then  fully  determined  by  F and  N(0). 

If  the  shock  is  sufficiently  weak  at  x,  t , (5.  36)  shows  that  the  sign  of  the 
shock  strength  e = e(x,t)  is  determined^  by  the  dissipation  inequality 
' . (4.  19)?  This  in  turn  deternaines  the  sign  of  the  jumps  across  the  shock 
in  other  field  quantities  such  as  the  density  (see  (5.  8)  ), 

If  the  compressible  material  at  hand  is  isotropic,  it  is  possible  to 
use  (5.  37)  (5.  22),  (5.  4)  and  (5.  19)  to  show  that  D assumes  the  following 
special  form: 

° = -ttJ'wjjj  - iKj  + t 

- JWiijC^  +|-WiuC^  . (5.38) 


If  N(0)  satisfies  (5.  20),  so  does  -N(0)  (see  (5.  17)  ).  Thus  if  x is  a non- 
contact  point  of  (as  is  assumed  here),  then  a vector  satisfying  (5.  20) 

may  be  assumed  to  satisfy  X.(5*  without  loss  of  generality. 

2 

Except  in  the  special  case  when  D = 0;  this  eventuality  will  not  be  considered 
here.  i j 

^See  [3]  for  details  in  special  cases.  | I 
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where 

^ . (5.  39) 

and  the  derivatives  of  W appearing  in  (5.  38)  are  to  be  evaluated  at  the 

+ + + 
values  I and  J appropriate  to 

6.  Incompressible  materials:  equilibrium  shocks  in  anti-olane  strain. 
A simpler  class  of  deformations  for  the  cylindrical  body  of  the 
preceding  section  is  that  of  anti-plane  strain,  in  which  the  displacement  vector 
satisfies 

^1  = U2  = 0 , u^=u(^Xj,X2^  . (6.1) 

A particle  is  therefore  displaced  only  axially,  and  the  amount  of  displacement 
depends  only  upon  the  position  of  the  particle  in  its  cross-section.  It 
follows  from  the  first  of  (2.  3)  that  the  deformation  gradient  tensor  F asso- 
ciated with  the  displacement  field  (6.  1)  has  components 

^Qrp  “ ' ^03  ~ ® ' ^30  “ a » ^33  = ^ • (^*  ^) 

From  (6.  2)  and  (2.  13)  one  obtains 

h " ^2  " ^ * J=detF  = l , (6.3) 

the  last  of  which  shows  that  anti-plane  strain  is  locally  volume-preserving. 

For  anti-plane  strain  the  field  equations  of  Section  2 can  be  reduced 
for  a compressible  material  to  three  differential  equations  for  the  single 
unknown  u , and  in  the  incompressible  case  to  three  differential  equations 


for  the  two  unknowns  u and  p.  It  is  thus  not  surprising  that  not  all  materials 
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— not  even  all  homogeneous,  isotropic  ones  — can  sustain  general  states  of 
anti-plane  strain  in  the  absence  of  body  forces?  One  class  of  incompressible 
materials  which  do  have  this  property  consists  of  those  whose  strain  energy 
densities  depend  only  on 

W = w(lj  , W(3)  = 0 . (6,  4)^ 

Such  materials  are  isotropic  and  include  the  neo-Hookean  material,  in  which 
W is  linear  in  , as  a special  case.  If  (6.  4)  holds,  the  governing  equations 
(2. 6),  (2.  19),  (6.  2),  (6.  3)  and  the  third  of  (2.  3)  can  be  consistently  reduced 
to  the  following  single  differential  equation^  for  the  out-of-plane  displacement 
u^Xi,X2)  : 


= 0 


on  ft  , 
o 


(6.  5) 


where,  as  in  the  preceding  section,  ft^  denotes  an  open  cross-section  of 
the  undeformed  cylinder,  and  the  prime  indicates  differentiation  with  respect 
to  the  argument  of  W. 


A necessary  and  sufficient  condition  that  a homogeneous,  isotropic,  in- 
compressible elastic  material  admit  general  states  of  anti-plane  strain  is 
derived  in  [18];  the  corresponding  condition  for  compressible  materials  is 
considerably  more  restrictive  and  may  be  found  in  [i9j. 

^The  requirement  that  W should  vanish  in  the  undeformed  state  is  a con- 
venient normalization. 

^Details  of  this  reduction  are  included  in  [18];  it  is  assumed  that  the  tractions 
on  the  lateral  surface  of  the  cylinder  are  independent  of  x^.  See  also  [6]. 
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The  components  of  nominal  stress  £ and  true  stress  ' are  given  by 

- "03  =°3<.  ' 

^0(3  = = °33  = ° ' (6.6 

"33  = ^ • 


If  u = kx^,  where  k is  constant,  then  (6,  5)  is  satisfied  and  the 

2 

deformation  is  one  of  simple  shear.  The  relation  between  the  shear  stress 
■'■^3  “ and  the  amount  of  shear  k is  found  from  (6.  6)  to  be 


T = 2w'(^3  + k^)k  (6.  7) 

in  simple  shear;  (-1  = 2W  (3)  is  the  shear  modulus  for  infinitesimal  defor- 
mations. The  shear  modulus  at  arbitrary  shear  k is  2W^^3  + k^ j and 
is  assumed  to  be  positive  for  all  k. 

The  differential  equation  (6.  5)  is  found  to  be  elliptic  at  a solution 
u and  at  a point  ^x^,  x^  ) if 


(6.8) 


From  (6.  7),  this  is  equivalent  to  the  condition  that  the  curve  of  shear  stress 


The  hydrostatic  pressure  p occurring  in  (2.  19)  has  been  eliminated  in 
terms  of  |Vu|  in  the  course  of  reducing  the  field  equations  to  (6.  5).  See 
[ 18]  or  [6]. 

2 

Simple  shear  can,  of  course,  be  sustained  without  body  forces  in  any 
elastic  material.  * 
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versus  amount  of  shear  shall  have  positive  slope^  at  k = j^u^x^.x^yl. 

For  anti-plane  strain,  the  components  of  the  energy-momentum  tensor 
are  found  from  (4.  11),  (6.2)  and  (6.6)  to  be 


(6.9) 


When  taken  over  closed  cylindrical  surfaces  T of  unit  length  which  are 
coaxial  with  K,  the  integrals  of  (4.  13)  reduce  to 


= - ^ Zw'u.pNpds  , 


(6.  10) 


(6.  11) 


where  is  the  closed  curve  forming  the  boundary  of  the  cross-section  of 
r,  and  s is  arc -length  on  If  encloses  only  points  of  where 

u is  twice  continuously  differentiable , one  has  the  conservation  laws 

= 0 , J3  = 0 . (6.  12) 

as  can  be  verified  directly  with  the  aid  of  (6.  5)  and  the  divergence  theorem. 
The  function  H of  (4.  17)  reduces  because  of  (6.9)  to 

H = - W(3  + |7u|2}  + 2W'(3  + IVul^j  (1^/  . (6.  13) 


^For  the  neo-Hookean  material,  W ^3  + = U = constant,  (6.  5)  is  L.aplace's 

equation,  and  (6.  8)  holds  for  all  k. 

W and  are  to  be  evaluated  at  3 + l_^u|  . 
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where  9u/9N 

A special  subclass  of  hypothetical  incompressible  isotropic  materials^ 
which  satisfy  (6.  4)  are  those  for  which 

=roj['  *^'1  ‘}  ■ 

where  l-i>0,  b>0  and  n>0  are  material  constants.  The  T _ relation 

for  this  class  of  materials  is 

T = la  ^1  + . (6.  15) 

Graphs  of  T vs  k for  various  values  of  n are  shown  in  Fig.  2.  The  el- 
lipticity  condition  (6.8)  is  always  satisfied  if  b>0,  n^l/2.  If  b>0,0<n<l/2, 

ellipticity  of  (6.  5)  fails  if 


l~u(x,.X2)l 

It  will  be  assumed  henceforth  that  0<n<l/2.  For  the  special  materials  char- 
acterized by  (6.14),  formula  (6.  13)  for  H reduces  to 


(6.  17) 


For  weak  solutions  in  anti-plane  strain  which  involve  equilibrium 
shocks,  both  the  material  and  the  spatial  shock  surfaces  are  cylindrical,  and 


These  materials  have  been  considered  in  connection  with  a crack  problem 
in  [ 6]. 

2 

One  always  has  in  locally  volume -preserving  deformations. 
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a material  shock  intersects  the  cross-section  ft  in  a curve  C.  The  traction 

o 

jump  condition  (3.  3)  can  be  seen  from  (6.  6)  tc  reduce  to 


[w'(.3  < Tn^  = 


on  C , 


(6.  18) 


while  the  conditions  (3.  5)  for  continuity  of  the  tangential  displacement  gradients 


become 


(“•ol!  <-0  --  0 


on  C , 


(6.  19) 


where  N and  L are  unit  vectors  normal  and  tangent  to  C,  with  N pointing 
into  the  positive  side  of  C. 

From  (6.  19),  one  concludes  that  there  is  a scalar  g such  that 


u,  ^ = u,  ^ + gN^ 


on  C 


(6.  20) 


Let  £ be  introduced  by 


1 - e = 


W'  {3  + l^ul^j 
W'(3  + 


(6.21) 


thus  € is  the  relative  change  in  the  shear  modulus  across  the  shock.  Since 
it  has  been  assumed  that  W >0,  one  has  6<  1.  The  traction  continuity  con- 
dition (6.  18)  can  be  written  as 


(1  on  C 


(6.  22) 


whence,  from  (6.  20) 


Equation  (6.  20)  is  the  analog  in  anti-plane  strain  of  (5.  10)  for  plane  strain. 
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g 


From  (6,  20),  (6.  23)  one  has 


(6.  23) 


on  C , 


(6.24) 


so  that  (6.  21)  reduces  to 


w'(3  + k ) = (1  - e)w'{3  + k - 


8u  ^2 

aN/ 


(2  - e)ej. 


on 


(6.25) 


where  the  abbreviations 

k = (6.26) 

have  been  introduced.  If  ^u  , and  therefore  Ic  , is  regarded  as  given, 

(6.  25)  furnishes  one  equation  for  the  determination  of  du/ dN  = * N as  a 

function  of  €.  According  to  (6.  20)  and  (6.  23),  one  would  thus  expect  to 
find  a one-p>arameter  family  of  gradients  _|Ju  which  can  be  "connected" 

to  ^u  via  the  shock  conditions^  e — the  shock  strength  — is  to  be  chosen  as 

2 

the  indexing  parameter. 

For  weak  shocks  in  anti-plane  strain,  Ic  1,  and  an  analysis 
parallel  to,  but  much  simpler  than,  that  carried  out  in  the  preceding 
section  shows  that  the  limiting  value  as  e "♦  0 of  du/ dN  satisfies 

^It  is  possible  to  show  that  the  local  existence  of  a solution  du/dN  of  (6.  25) 
for  some  e ^ 0 requires  the  material  to  be  such  that  the  ellipticity  condition 
(6.  8)  should  fail  for  some  value  of  the  amount  of  shear  k. 

^Note  that  the  physical  meaning  of  e in  anti-plane  strain  is  not  the  same  as 
that  of  its  counterpart  for  plane  deformations;  compare  (6.  21 ),  (5,  8). 
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/Ou  f -W'(3  + k^) 

, 0 ■ 2w"(i  t k-") 


(6.  27) 


This  result,  in  turn,  can  be  used  to  show  that  the  limiting  direction  of  the 
shock  as  €">0  is  that  of  a characteristic  curve  of  the  partial  differential 
equation  (6.  5).  There  are  in  general  two  characteristic  directions  associated 
with  a point  (x^.x^)  at  which  ^u(x^,X2)  is  such  that  ellipticity  has  failed,  and 
these  characteristic  directions  are  symmetric  with  respect  to  ^u(Xj,X2). 

By  carrying  out  the  appropriate  small  - je  [ approximations  using 
(6.  13)  and  (6.  25),  one  obtains  the  following  weak-shock  formula  for  [H]j 


m!  Km)e  =0 


, 3+  \ 2 


.W 


4 O(e^)  as  e -*0 


(6,  28) 


where  W*  and  w"'  are  to  be  evaluated  at  3 + k , and  k is  given  in 
(6.  26).  In  this  form,  the  approximation  (6.  28)  is  analogous  to  the  plane 
strain  result  <5.  36),  (5.  38).  In  the  present  circumstances,  (6.  27)  can  be 
used  to  further  reduce  (6.  28)  to 


/^u^^ 


3(w")^ 


/ //< 

- w w 


e + 0(e  ) as  e -♦0 


(6.29) 


where,  again,  the  derivatives  of  W are  evaluated  at  3 + l^u|^. 

Now  consider  a quasi-static  time -dependent  family  of  equilibrium 
shocks  in  anti-plane  strain,  and  let  the  shock  strength  at  time  t associated 
with  this  family  be  e(t).  If  the  shock  is  sufficiently  weak  at  time  t,  and  if 
the  dissip>ation  inequality  (4.  19)  holds,  then  (6.29)  shows  that 
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3[w"(lj)]^  - w'(lj)W"'(lj)  e^(t)  so  , 


(6.  30) 


where  = 3 + k (t). 


For  the  special  material  characterized  by  the  strain  energy  density 
(6.  14),  one  verifier  easily  that 


3(W")‘^  - w'w'" 


' 2n) 

4ir 


1 + - (I,  - 3) 
n 1 


2n-4 


(6. 31) 


Since  b>0,  0-  n<l/2,  this  quantity  is  negative,  and  (6.30)  thus  requires  that 

G(t),  which  is  always  less  than  unity,  shall  also  be  nonnegative  for  this 

material,  at  least  for  sufficiently  weak  shocks.  Reference  to  (6.  21)  then 

/ 2 

shows  that  the  shear  modulus  2W  (3  + k ) is  smaller  on  the  positive  — 
or  "upstream"  — side  of  the  shock  than  on  the  downstream  side.  Moreover, 

(6.24)  shows  that,  for  a dissipative  weak  shock  in  the  particular  material 
under  consideration,  the  downstream  side  carries  the  smaller  value  of 
|^u|^,  since  0<e<l. 

For  the  special  strain  energy  density  (6.  14),  the  basic  equation 

(6.25)  governing  the  shock  transition  can  in  fact  be  solved  explicitly  and 

the  conclusions  above  can  be  valida.ted  for  all  equilibrium  shocks  — not 

merely  weak  ones.  The  result  of  solving  (6.  25)  in  this  case  is 

1 


n / €(2  - €) 


1 - n 


e < 1 . 


(6.  32) 


Recall  that  the  velocity  V of  the  quasi-statically  moving  shock  C points 
into  the  "positive"  side  of  C. 


one  must  have 


1 


+ 2 


(6.  34) 


if  a shock  of  strength  e is  to  exist  in  the  presence  of  a displacement  grad- 
+ 

ient  of  magnitude  k on  the  positive  side  of  the  shock.  Thus,  for  a given  value 

+ 

of  the  material  constants  b and  n,  only  those  points  in  the  e,  k — plane 
which  lie  above  the  curve  in  Fig.  3 — the  admissible  region  — correspond  to 
equilibrium  shocks  in  anti-plane  strain. 

Finally,  one  can  compute  [H]^  exactly  for  the  special  material  at 
hand  from  (6.  13),  (6.  14)  and  (6.  32).  One  finds 


(6.35) 


where 


=FTT  l^n  - 2 + € + (1  - c)^  - " + (1  - 2n){l  - e)‘  " " |.  (6.36) 

It  can  then  be  shown  that  f(e)<0  for  0<e<l,  f(0)  =0,  and  f(e)>0  for 
- CD  <€  <0.  The  dissipation  inequality  (4.  19)  would  then  again  require 
0<€  < 1 for  all  (not  just  weak)  shucks.  That  portion  of  the  admissible 
region  of  Fig.  3 which  corresponds  to  dissipative  shocks  is  indicated  by 
shading.  Note  that  ellipticity  is  always  lost  on  the  upstream  side  of  the  shock 
when  the  dissipation  condition  holds. 
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1 


7.  A thermodynamic  basis  for  the  disai^ati^n  inequality. 

The  present  section  is  concerned  with  a thermodynamic  argument^ 
in  support  of  the  dissipation  inequality  (4.  19).  For  simplicity,  only  com- 
pressible materials  will  be  considered. 

Consider  a material  characterized  by  an  internal  energy  e = e(F,  T]) 
(measured  per  unit  of  mass)  which  depends  only  on  the  deformation  gradient 
tensor  F = F(x,  t)  and  the  entropy  per  unit  mass  11  = i1(x,  t),  and  for  which 

the  nominal  stress  tensor  ^=£(x,t)  and  the  temperature  9 =0(x,  t)>0  are 

.2 

given  by 


0.  . 
IJ 


= P 

*^o  9F.. 


(7.  1) 


Here  is  the  (constant)  mass  per  unit  undeformed  volume,  and,  as  in 

the  preceding  sections,  x is  the  position  vector  to  a particle  in  the  unde- 
formed state.  The  internal  energy  density  e is  assumed  to  be  infinitely 
differentiable  with  respect  to  F and  ti,  and  also  to  be  such  that  the 
second  of  (7.  1)  is  uniquely  invertible  to  give 

= ^(F,  0)  (7.  2) 


as  an  infinitely  differentiable  function  of  F and  0 for  every  nonsingular 
F and  every  0 > 0,  Let 

♦ = ♦(F,0)  = e(F,f|(F,0))  - 0f|(F,0)  (7.  3) 


I am  indebted  to  Professor  James  R.  Rice  for  the  substance  of  the  argument 
given  here. 

In  the  terminology  of  580  of  [ 20],  such  a material  is  called  perfect. 
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be  the  free  energy  per  unit  mass.  One  shows  easily  from  (7.  I)  - (7.  3) 
that 


“ij = ■ '’•“I 

u 

As  in  section  3,  it  is  assumed  that,  for  each  time  t in  the  interval 
of  interest,  the  displacement  field  u is  continuous  and  piecewise  continuously 
differentiable  on  the  region  tt  occupied  by  the  undeformed  body,  so  that 
F = ^u  is  piecewise  continuous  on  R for  each  t.  It  is  also  assumed  that 
Tl  is  piecewise  continuous  on  R for  each  t.  Finally,  ja,  F and  T]  are 
presumed  to  be  continuously  differentiable  in  t for  each  x in  R.  It 
follows  from  (7.  1)  that  the  smoothness  properties  of  0 and  0 are  the 
same  as  those  of  F and  t) . 

Let  & be  an  arbitrary  regular  subregion  of  R,  and  let  be  the 
image  at  time  t of  £ under  the  motion  ^ ^ + iu(x,t).  The  balance  of 

energy  is  postulated  as  follows: 

^ J Po®^^  ^ I J =j£-vdA+  Jh-NdA  , (7.5) 

i £ £ i d£  8£ 

for  every  £.  Here  v = v(x,  t)  is  the  particle  velocity,  ja  =ON  is  the 

sjc 

nominal  traction  acting  on  the  boundary  3£^  of  £^  , and  h is  the  heat 
flux  (measured  per  unit  area  of  d£)  acting  on  3£^  . Thus  the  left  side 
represents  the  rate  of  increase  of  the  sum  of  internal  and  kinetic  energy, 

while  the  right  side  is  the  sum  of  the  rate  of  work  of  the  tractions  and  the 

% 1 

rate  of  heat  flow  into  £^  across  its  boundary. 

^Body  forces  and  internal  heat  sources  are  assumed  to  be  absent. 
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The  rate  of  entropy  production  T (t)  associated  with  the  motion  is 
defined  by 

J - J^h-NdA  . (7.6) 

& d& 

It  is  further  postulated  that 

r(t)2:0  (7.7) 

for  every  &,  for  every  motion  and  for  each  time  t. 

Suppose  now  that  the  temperature  is  uniform  throughout  the  body 
and  independent  of  time  0(x,  t)  = 9^  = constant  for  every  x and  W and  all 
t.  By  eliminating  the  heat  flux  integral  between  (7.  5)  and  (7.  6)  and  making 
use  of  the  definition  (7.  3),  one  finds  that 


Thus  the  entropy  production  rate  is  proportional  to  the  excess  of  the  rate  of 
work  of  the  forces  external  to  over  the  rate  of  increase  of  the  sum  of  the 
free  energy  and  the  kinetic  energy. 

A comparison  of  the  first  of  (7.  4)  with  (2,  9)  shows  that  the  material 
presently  under  consideration  may  be  identified  with  the  (compressible) 
elastic  material  of  section  2 by  writing 

= Pq'*' ®o^  ' 

as  long  as  the  motion  is  isothermal.  For  isothermal  quasi-static  motions 

of  the  kind  considered  in  section  4,  one  neglects  the  contribution  of  the  ! 

i' 

t 

I 
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kinetic  energy  in  (7.  8)  and  thus  obtains,  with  the  aid  of  (7.  9),  the  result 


o 


1 s • vdA  - U 

I J ~ ~ 


(7.  10) 


where 


U(t)  = J W{F(x.t))dV  . (7.  11) 

£ 

Now  suppose  that  the  isothermal,  quasi-static  motion  involves  a 
shock  with  the  properties  described  in  section  4.  If  the  material  shock 
surface  intersects  £,  the  (4.  16)  applies  and  can  be  used  to  write  (7.  10) 
in  the  form 


(7.  12) 


where  is  the  part  of  lying  in  £.  The  postulate  (7.7)  of  nonnegative 
entropy  production  rate  then  leads  immediately  to  (4.  18),  from  which  the 
dissipation  inequality  (4.  19)  follows. 

If  no  shock  is  present,  then  [H]  = 0 and  (7.  12)  gives  r(t)  = 0.  Thus, 
although  smooth  quasi -static  motions  of  an  elastic  body  do  not  result  in 
dissipation,  weak  solutions  involving  equilibrium  shocks  generally  do. 
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FIGURE  3.  ADMISSIBLE  REGION  FOR  EQUILIBRIUM  SHOCKS  IN 
ANTI-PLANE  STRAIN  FOR  AN  INCOMPRESSIBLE 
MATERIAL  CHARACTERIZED  BY  (6.14)  WITH  0<n<l/2. 
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